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Abstract We discuss the ground state properties of matter in outer and inner crusts 
of neutron stars under the influence of strong magnetic fields. In particular, we 
demonstrate the effects of Landau quantization of electrons on compositions of neu- 
tron star crusts. First we revisit the sequence of nuclei and the equation of state 
of the outer crust adopting the Baym, Pethick and Sutherland (BPS) model in the 
presence of strong magnetic fields and most recent versions of the theoretical and 
experimental nuclear mass tables. Next we deal with nuclei in the inner crust. Nu- 
clei which are arranged in a lattice, are immersed in a nucleonic gas as well as a 
uniform background of electrons in the inner crust. The Wigner-Seitz approxima- 
tion is adopted in this calculation and each lattice volume is replaced by a spherical 
cell. The coexistence of two phases of nuclear matter - liquid and gas, is consid- 
ered in this case. We obtain the equilibrium nucleus corresponding to each baryon 
density by minimizing the free energy of the cell. We perform this calculation using 
Skyrme nucleon-nucleon interaction with different parameter sets. We find nuclei 
with larger mass and charge numbers in the inner crust in the presence of strong 
magnetic fields than those of the zero field case for all nucleon-nucleon interactions 
considered here. However, SLy4 interaction has dramatic effects on the proton frac- 
tion as well as masses and charges of nuclei. This may be attributed to the behaviour 
of symmetry energy with density in the sub-saturation density regime. Further we 
discuss the implications of our results to shear mode oscillations of magnetars. 
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1 Introduction 



Magnetised Neutron star crusts 



Neutron star crust is a possible site where neutron rich heavy nuclei might reside. 
Extreme physical conditions exist at the crust of a neutron star. The temperature is 
^ lO"^ K and the density varies from 10^ - lO''* g/cm^ there. Recently, it was ob- 
served that certain neutron stars called magnetars had surface magnetic fields ^ lO'^ 
G. The internal fields could be several times higher than the surface fields of mag- 
netars. Soft gamma repeaters (SGRs) are suitable candidates for magnetars ||T]|2]|3|. 
Giant flares were observed from SGRs in several cases. Those giant flare events are 
thought to be the results of star quakes in magnetars. This might be attributed to 
the tremendous magnetic stress due to the evolving magnetic field leading to cracks 
in the crust. Quasi-periodic oscillations discovered in three giant flares are the evi- 
dences of torsional shear mode oscillations in magnetar crusts. 

Such strong magnetic fields in magnetars are expected to influence charged parti- 
cles such as electrons in the crust through Landau quantization. The effects of strong 
magnetic fields on dense matter in neutron star interior were studied earlier |4jl5j. 
It was also noted that atoms, molecules became more bound in a magnetic field |0 . 
In this article, we discuss the effects of strongly quantising magnetic fields on com- 
positions and equation of state of the ground state matter in neutron star crusts and 
its connection to torsional shear mode oscillations. 

We organise the article in the following way. Neutron star crusts in strong mag- 
netic fields are described in sections 2, 3 and 4. Torsional shear mode oscillations of 
magnetars are discussed in section 5. Finally, we summarise in section 6. 



2 Crusts in Strong Magnetic Fields 



We investigate compositions and equations of state (EoS) of outer and inner crusts 
in strong magnetic fields. Nucleons are bound in nuclei in the outer crust. Nuclei 
are immersed in a uniform background of electron gas which becomes relativistic 
beyond 10^ g/cm^. Neutrons start to drip out of nuclei at higher densities. This is the 
beginning of the inner crust. In this case, nuclei are embedded both in electron and 
neutron gases. Magnetic fields may influence the ground state properties of crusts 
either through magnetic field and nuclear magnetic moment interaction or through 
Landau quantisation of electrons. In a magnetic field ^ 10^^ G, magnetic field and 
nuclear magnetic moment interaction would not produce any significant change. 
However such a strong magnetic field is expected to influence charged particles 
such as electrons in the crust through Landau quantization. Our main focus is to 
study the effects of Landau quantisation on the ground state properties of neutron 
star crusts. Later we discuss shear mode frequencies using our results of magnetised 
neutron star crusts. 
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2.1 Landau Quantisation of Electrons 



We consider electrons are noninteracting and placed under strongly quantising mag- 
netic fields. In the presence of a magnetic field, the motion of electrons is quantized 
in the plane perpendicular to the field. We do not consider Landau quantisation of 
protons because magnetic fields in question in this calculation are below the critical 
field for protons. However, protons in nuclei would be influenced by a magnetic 
field through the charge neutrality condition. We take the magnetic field (^) along 
Z-direction and assume that it is uniform throughout the inner crust. If the field 
strength exceeds a critical value Be = m'i/e ~ 4.414 x lO'^G, then electrons be- 
come relativistic ||6|. The energy eigenvalue of relativistic electrons in a quantizing 
magnetic field is given by 

Ee{v,Pz) = [pl + ml + 2eBv]''^ , (1) 

where is the Z-component of momentum, v is the Landau quantum number. 
The Fermi momentum of electrons, pPev' is obtained from the electron chemical 
potential in a magnetic field 

PF,,,= [^le^-m;-2eBvf\ (2) 
The number density of electrons in a magnetic field is calculated as 

Vmm- 

= ^ L SvPF,,y , (3) 

v=0 

where the spin degeneracy is gy = I for the lowest Landau level (v = 0) and gv ^2 
for all other levels. 

The maximum Landau quantum number (Vmm) is obtained from 



2 2 

V,na. - 

The energy density of electrons is. 



(4) 



eB / , , , PF,y+^J■e \ 

4^ v=o V ^{ml+2eBv)) 
Similarly the pressure of the electron gas is determined by 

= £ PF^.ylle - {ml + 2eBv)\n ^ I" . (6) 

47r2 \ ^{ml+2eBv) / 
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Magnetised Neutron star crusts 




Fig. 1 Gibbs free energy per nucleon is plotted with mass density for zero magnetic field (left 
panel) and — (right panel). Equilibrium nuclei are shown with solid symbols in both panels. 



3 Magnetic BPS Model of Outer Crust Revisited 

We describe the BPS model in the presence of strong magnetic fields B ~ lO'^G 
to determine the sequence of equilibrium nuclei and the equation state of the outer 
crust Q |8). Nuclei are arranged in a bcc lattice in the outer crust. Here we adopt 
the Wigner-Seitz (WS) approximation and replace each lattice volume by a spher- 
ical cell which contains one nucleus at the center Further each cell is to be charge 
neutral such that equal numbers of protons and electrons are present there. The 
Coulomb interaction among cells is neglected. An equilibrium nucleus (A, Z) at a 
given pressure P is obtained by minimising the Gibbs free energy per nucleon with 
respect to A and Z. In this calculation, we modify the magnetic BPS model includ- 
ing the finite size effect in the lattice energy and adopting recent experimental and 
theoretical mass tables. The total energy density of the system is given by 

E,„,=nN{WN+WL) + ee. (7) 

The energy of the nucleus (including rest mass energy of nucleons) is 

Wn = m„ (A - Z) + mpZ - M , (8) 

where is the number density of nuclei, b is the binding energy per nucleon. 
Experimental nuclear masses are obtained from the atomic mass table compiled 
by Audi, Wapstra and Thibault |j9l- For the rest of nuclei we use the theoretical 
extrapolation of MoUer et al. ifTOl . Wl is the lattice energy of the cell and is given by 

Here rc is the cell radius and ~ r^A^I^ (ro ~1.16 fm) is the nuclear radius. The 
first term in Wl is the lattice energy for point nuclei and the second term is the 
correction due to the finite size of the nucleus (assuming a uniform proton charge 
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distribution in the nucleus). Further £e is the electron energy density as given by 
Eq.(|5]l and P is the total pressure of the system given by 



where Pe is the pressure of electron gas in a magnetic field as given by Eq.®. 
The Gibbs free energy per nucleon is 



where n is the total baryon number density. 

At a fixed pressure P, we minimise g varying A and Z of a nucleus. The sequence 
of equilibrium nuclei and their corresponding free energies are shown in FiglT] Here 
we define = B/Bc- The left panel shows results for B ~Q and the right panel 
corresponds to = 10^. It is evident from the figure that some nuclei disappear and 
new nuclei appear under the influence of strong magnetic fields. It is attributed to the 
phase space modification of electrons due to Landau quantisation which enhances 
the electron number density JS) . 



4 Inner Crust in Quantizing Magnetic Fields 

Now we describe the ground properties of matter of inner crusts in presence of 
strong magnetic fields using the Thomas-Fermi (TF) model at zero temperature. In- 
ner crust nuclei are immersed in a nucleonic gas as well as a uniform background 
of electrons. Furthermore, nuclei are arranged in a bcc lattice. As in the case of 
outer crust, we again adopt the Wigner-Seitz (WS) approximation in this calcu- 
lation. Here each cell is taken to be charge neutral and the Coulomb interaction 
between cells is neglected. Electrons are uniformly distributed within a cell. The 
system is in jS -equilibrium. We assume that the system is placed in a uniform mag- 
netic field. Though electrons are directly affected by strongly quantizing magnetic 
fields, protons in the cell are influenced through the charge neutrality condition ifTTl . 
The interaction of nuclear magnetic moment with the field is not considered because 
it is negligible in a magnetic field below lO'** G ||T2| . 

The spherical cell in the WS approximation does not define a nucleus. We exploit 
the prescription of Bonche, Levit and Vautherin |1T3] [141 to subtract the gas part 
from the cell and obtain the nucleus. It was shown that the TF formalism at finite 
temperature generated two solutions llTSi - one for the nucleus plus neutron gas and 
the other representing the neutron gas. The nucleus is obtained as the difference of 
two solutions. This formalism is adopted in our calculation at zero temperature as 
described below. 

The thermodynamic potentials for nucleus plus gas (NG) and only gas (G) phases 
are defined as l[T3l[T4l 




(10) 



g 



Eror+P _ WN+A/2>WL + Zlle 
n A 



(11) 



q=n.,p 
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where ji^ and are the free energy density, baryon chemical potential and 
number density, respectively. The nucleus plus gas solution coincides with the gas 
solution at large distance i.e. £2ng = ^G- Ths free energy which is a function of 
baryon number density and proton fraction (Yp), is defined as ifTTl 



(13) 



^{rigjp) = [■Jif + £c + £e\dr . 



The nuclear energy density is calculated using the Skyrme nucleon-nucleon in- 
teraction and it is given by lfT6l[T7l[T8l 
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and the effective nucleon mass 
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(15) 



where total baryon density is n = n„ + iip. 

The direct parts of Coulomb energy densities for the nucleus plus gas and gas 
phases follow from lfm[T9l 



— («^(r')-«.K 



(16) 



where and «p are proton densities in two respective phases. The exchange parts 
of coulomb energy densities are small and neglected in this calculation. 

The average electron chemical potential in a magnetic field given by Eq.© is 
modified to HIl 
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^e^[pF,Avf + ml + 2eBvf^-<V'ir)> , (17) 

where <V'^{r) > denotes the average single particle Coulomb potential and for both 
phases it is given by 

y'i'-)=J KV)-".]^^^^'. (18) 

The density profiles of neutrons and protons with or without magnetic fields are 
obtained by minimising the thermodynamic potential in the TF approximation 



dQ 



NG 



5Q., 



G 



5«G 



= 0, 

= 0, (19) 



with the condition of number conservation of each species from 

Zceii = I nf{r)dr , 

NG 



Nceii = J <''{r)dr , (20) 

where N^eii and Z^eii are neutron and proton numbers in the cell, respectively. 

We obtain the mass number A=N + Z and atomic number using the subtraction 
procedure as 



Z = j [n^,^{r)-n^^{r)\dv, 

f [n:!^ir)~n^{r)]dr. (21) 



Here we again obtain the equilibrium nucleus at each density by minimising the 
free energy of the nuclear cluster in the cell along with charge neutrality and j3- 
equilibrium conditions ifTTl . In the left panel of Fig|2] proton fraction is shown for 
B ~0 and = 10^. Protons are influenced by the Landau quantisation of electrons 
through charge neutrality condition. At lower densities, only the zeroth Landau level 
is populated by electrons whereas a few Landau levels are populated above density 
0.005 fm"^ for = 10"^ i.e. 4.414x lO'^ G. This is reflected in the proton fraction 
which rises hugely at lower densities and approaches to the zero field case at higher 
densities. Further we estimate the effects of different parameter sets of Skyrme in- 
teraction on the proton fraction. It is noted that the SLy4 set |[20l results in higher 
proton fraction due to the stiffer density dependence of the symmetry energy at sub- 
saturation densities than that of the SkM set. 

We exhibit mass and atomic numbers of equilibrium nuclei after subtraction of 
free neutrons as a function of average baryon density in the right panel of Fig|2l 
Results are obtained for B = and = 10^. Besides SkM and SLy4 parameter 
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Fig. 2 Proton fraction (left panel) and mass and atomic numbers of equilibrium nuclei (right panel) 
are plotted with average baiyon density for different magnetic field strengths and Skyrme interac- 
tion parameter sets. 



sets, we also exploit Sk272 II2TI parameter set for this calculation. In all three cases, 
mass and atomic numbers are higher than zero field cases as long as only the zeroth 
Landau level is populated. However, the situation is changed at higher densities 
when electrons jump from the zeroth Landau level to the first level. This leads to 
jumps in mass and atomic numbers in nuclei as noted for the SLy4 set. Further, 
the variation of parameters for nucleon-nucleon interaction affects mass and atomic 
numbers of nuclei as it is evident from the figure. We also note that the free energy 
of the ground state matter in strong magnetic fields is reduced and becomes more 
bound compared with the field free case. 



5 Shear Mode Oscillations in Magnetars 

Giant x-ray flares caused by the tremendous magnetic stress on the crust of magne- 
tars were observed in several cases. Star quakes associated with these giant flares 
excite seismic oscillations. Quasi-periodic oscillations (QPOs) were found in the 
decaying tail of giant flares from SGR 1900H-14 and SGR 1806-20. Those QPOs 
were identified as shear mode oscillations of magnetar crusts |l3]|22l. Frequencies 
of the observed QPOs ranged from 18 Hz to 1800 Hz. 

Shear mode frequencies are sensitive to the shear modulus of neutron star crust. 
The shear modulus is again strongly dependent on the composition of neutron star 
crust. It might be possible to constrain the properties of neutron star crusts by study- 
ing the observed frequencies of QPOs. Torsional shear mode oscillations were in- 
vestigated both in Newtonian gravity Il23ll24l and general relativity 125, 26, , 271 . In 
both cases, it was assumed that the magnetised crust was decoupled from the fluid 
core. 

Here we describe the calculation of shear mode frequencies adopting the model 
of Sotani et al. Il26l . In this case, we study torsional shear oscillations of spherical 
and non-rotating relativistic stellar models. The metric used here has the form. 
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ds^ = -e^'^dt^ + e^'^dr + r^{d9- + sin^ed(j)^) . (22) 

The equilibrium models are obtained by solving Tolman-Oppenheimer-Volkoff 
equation. Next the equilibrium star is assumed to be endowed with a strong dipole 
magnetic field 1261 . The deformation in the equilibrium star for magnetic fields 
is neglected. Torsional shear modes are the results of material velocity 
oscillations. These modes are incompressible and do not result in density perturba- 
tion in equilibrium stars. Consequently, this leads to negligible metric perturbations 
and justifies the use of the relativistic Cowling approximation 1261 . The relevant 
perturbed matter quantity for shear modes is the ^-component of the perturbed four 
velocity du'l' ^ 

dJ = e-1'd,iV(t,r)-^dgPi(cose) , (23) 
sinO 

where (9, and dg correspond to partial derivatives with respect to time and 0, respec- 
tively. Pi {cos ) is the Legendre polynomial of order £ and iV{t,r) is the angular dis- 
placement of the matter The perturbation equation is obtained from the linearised 
equation of motion. Finally, we estimate eigenfrequencies by solving two first order 
differential equations Eq.(69) and (70) of Sotani et al. Il26l . 

Now we study the dependence of shear mode frequencies on the compositions 
of magnetised crusts which are already described in sections 3 and 4. Earlier calcu- 
lations were performed with non-magnetised crusts [|3j |26] |27] |28l . One important 
input for the shear mode calculation is the knowledge of shear modulus of the mag- 
netised crust. Here we adopt the expression of shear modulus as given by 1291 [30l 

^ = 0.1194 ' , (24) 
a 

where a = 3/(4 7r«,), Z is the atomic number of a nucleus and is the ion density. 
This zero temperature form of the shear modulus was obtained by assuming a bcc 
lattice and performing directional averages llSTI . Later the dependence of the shear 
modulus on temperature was investigated with Monte Carlo sampling technique 
by Strohmayer et al. ll30l . However we use the zero temperature shear modulus of 
Eq.(l24]i in this calculation. 

We calculate the shear modulus using Eq.(l24]i and the compositions and equa- 
tions of state of magnetised crusts obtained in Sec. 3 and 4. This is shown as a 
function of normalised distance with respect to radius (R) of the star for different 
field strengths B ~ 0, = 10^ and = 10^ and a neutron star mass of 1.4 Mq in 
the left panel of Figl3] Shear modulus increases initially with decreasing distance 
and drops to zero at the crust-core boundary. For = 10"* i.e. 4.414x10'^ G or 
more, the shear modulus is enhanced appreciably compared with the zero field case. 

It was argued that shear mode frequencies are sensitive to shear modulus [3 28 1. 
We perform our calculation for shear mode frequencies using the model of Sotani et 
al. ||26| and the shear modulus of magnetised crusts as described above. We calculate 
fundamental shear mode frequencies for a neutron star mass of 1 .4 Mq as well as 



10 



Magnetised Neutron star crusts 




Fig. 3 Shear modulus is plotted as a function of noiTnalised distance for different magnetic field 
strengths (left panel) and shear mode frequencies are plotted with different £ values for a neutron 
star mass of 1.4 Mq and B = 4x 10''* G (right panel). 



magnetic fields as high as 4.414 x 10 G. When we compare those frequencies 
involving magnetised crust with those of the non-magnetised crust, we do not find 
any noticeable change between two cases. For SGR 1900+14 having Z? = 4 x lO'^ 
G and a neutron star mass of 1.4 Mq, we show in the right panel of Figl3]that 
the observed QPO frequencies match nicely with frequencies estimated using our 
magnetised crust model. Further we observe that the first radial overtones calculated 
with our magnetised crust model have higher frequencies than those calculated with 
the non-magnetised crust model. This is in agreement with the prediction that the 
radial overtones are susceptible to magnetic effects ll23l . 



6 Summary 

We have constructed the model of magnetised neutron star crusts and applied it 
to shear mode oscillations of magnetars. In particular, we highlighted the effects 
of strongly quantising magnetic fields on the properties of ground state matter of 
outer and inner crusts in this article. It is noted that compositions and equations 
of state of neutron star crusts are significantly altered in strong magnetic fields. 
Consequently, shear modulus of the crust which is sensitive to the compositions of 
crusts, is enhanced. We have observed that our model of the magnetised crust might 
explain the observed shear mode frequencies quite well. 
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